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Abstract. The paper investigates how correlations can completely specify a uniformly dis- 
crete point process. The setting is that of uniformly discrete point sets in real space for 
which the corresponding dynamical hull is ergodic. The first result is that all of the essential 
physical information in such a system is derivable from its n-point correlations, n — 2,3, ... . 
If the system is pure point diffractive an upper bound on the number of correlations required 
can be derived from the cycle structure of a graph formed from the dynamical and Bragg 
spectra. In particular, if the difi'raction has no extinctions, then the 2 and 3 point correlations 
contain all the relevant information. 



1. The Setting 

1.1. Quasicrystals and dynamical systems. The defining feature of physical cyrstals 
and quasicrystals is the prominent appearance of Bragg peaks in their diffraction diagrams. 
Mathematically the diffraction is a positive measure and the Bragg peaks comprise the pure 
point component of this measure. In a 'perfect' crystal or quasicrystal, the diffraction should 
be entirely pure point, and that is the situation that we shall assume here. We shall simply 
refer to these as quasicrystals in the sequel. 

The diffraction does not, on its own, determine the internal structure of the quasicrystal 
that created it. However, the diffraction does determine the 2-point correlation, which is its 
Fourier transform. Under appropriate conditions (assumptions A(i),(ii) below), knowledge of 
all the correlations (2-point, 3-point, etc.) does determine the internal structure (Thm. [L^ . 
The primary objective of this paper is to explore the details behind this in the case that 
the diffraction is a pure point measure. In this case, under fairly mild conditions, one does 
not at all need the entire set of correlations. In fact in the best situation, where there 
are no extinctions in the Bragg spectrum (this term is explained below), the 2- and 3-point 
correlations alone (in fact just the 3-point correlations) are enough to determine the structure 
(Thm.[2S]and its Corollary). 

Our approach here is to use a setting familiar from statistical mechanics and from the theory 
of tilings and long-range aperiodic order. Rather than deal with a single quasicrystal A, we 
work instead with translation invariant families of them. The intuition is that such a family, 
X, will consist of all those quasicrystals which are in some sense locally indistinguishable from 
one another, or which cannot be isolated from one another by the physical considerations at 
hand. As for the individual quasicrystals, we model these simply as uniformly discret^ll point 
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^A subset A of d-dimensional space R'' is uniformly discrete (or more specifically r-uniformly discrete) if 
for some r > and for all x,y £ A with x ^ y, \x — y\ > r. 
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sets in space, with the points representing the positions of the atoms. Since the dimension 
does not play any special role here, we work in general d-dimensional space M!^. 

Thus our main result deals with general ergodic uniformly discrete point processes. It 
is always the case that the dynamical spectrum is generated, as a group, by the diffraction 
spectrum. But here we prove that if it has pure point spectrum and the dynamical spectrum is 
expressible as a sum of finitely many copies of the diffraction spectrum then the point process 
is determined by its n-point-cor relations for some finite n. In the context of aperiodic order 
this becomes particularly relevant as it gives a way of assessing the degree of 'complexity' of 
the long range order. 

1.2. Background. It may be of interest to briefly discuss the reasons that such an elaborate 
formalism is relevant to what might seem a fairly straightforward exercise in spectral theory. 
Pure point diffraction from aperiodic structures was not predicted, either by mathematicians 
or crystallographers. When it did appear, both in tiling theory and experimentally in the 
discovery of aperiodic metallic alloys, the projection method was quickly utilized and it was 
generally believed that one could use standard techniques like the Poisson summation formula 
(applied to lattices in higher dimensions) to explain the diffraction. 

It was A. Hof who, in his much-cited papers [13^ I14j. showed that diffraction in aperiodic 
structures is not business as usual. The Bragg spectrum of an aperiodic material is not 
lattice-like as in the periodic case, but is typically dense in Fourier space. The problem is 
that for a countable aperiodic set A of scatterers in the Fourier transform i> of their Dirac 
comb u = "^xeA is not in general a measure. The sums involved diverge, even locally. This 
is in contrast to the lattice case. For this reason the theory of diffraction has developed by 
defining the diffraction as the Fourier transformation of the volume averaged autocorrelation. 
It is the 'quadratic' nature of autocorrelation which produces the necessary convergence, and 
this is now the standard approach to diffraction in the aperiodic case. 

Assume now that we are in the case of a countable and uniformly discrete set of scatterers. 
It was pointed out in [29] that a good way to study an aperiodic set was to follow ideas from 
statistical mechanics and form a compact space from its translation orbit, the so-called hull 
X. This is a dynamical system (with M*^ as the acting group) and allows one to use spectral 
theory. S. Dworkin [TT] then showed how the dynamical spectrum could be linked to the 
diffraction spectrum by using spectral measures. This linkage is now often called Dworkin's 
argument. 

It was the use of hulls and Dworkin's argument that first allowed rigorous proof of the pure 
pointedness of model sets (or cut and project sets). This is a good example of a situation where 
the result is seemingly clear from the Poisson summation formula, but on closer inspection 
one is confronted with divergent sums with no obvious mathematical meaning. 

The first proof by A. Hof, and its full generalization to all model sets by M. Schlottmann 
[3T] . of the pure pointedness of diffraction from model sets uses the fact that the hull is 
measure-theoretically a compact Abelian group, and so pure point, followed by Dworkin's ar- 
gument using spectral measures. By the way, unlike the periodic case, the hull is topologically 
not a group, and indeed its highly subtle topology has been the focus of many mathematicians 
recently (see \15\ [5] for reviews and further discussion). 

Dworkin's argument still left the precise connection between the diffraction and dynamical 
spectrum unresolved. Further developments on the hull and its connection to diffraction and 
to point processes were made in [5j and pj] where it is shown that under the assumption 
of ergodicity the autocorrelation of the points sets of X exist almost surely (in the sense 
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of the invariant measure /i on X) and almost surely are equal to the first moment of the 
Palm measure of X. In [8] this was extended to show that all the higher correlations of 
A G X exist almost surely and they completely determine /i. This fact is not generally true 
for point processes but here follows from the assumed uniform discreteness of the point sets 
under consideration. In the pure point case it comes pretty much for free from the spectral 
structure, as we see in the present paper. 

In [17j (see \\.2\ [2l [22] as well) it is proved that the diffraction is pure point if and only 
if the dynamical spectrum is pure point. This is remarkable since we know that in general 
the diffraction can fail to see great chunks of the dynamical spectrum - even the pure point 
part of the spectrum. In the diffraction/dynamics connection is made even more precise 
by showing that there is an isometric embedding of Lp'iW^^u}) into L^(X, where is 
the Fourier space with its diffraction measure lv. This embedding allows one to see how the 
eigenfunctions transfer across. Precisely, each Bragg peak located at position k gives rise to 
an eigenfunction whose value at the point set A € X is almost surely 



We mention this, first because the convergence of this sum away from is precisely what 
is called the Bombieri- Taylor conjecture in [131 114j . (which was precisely that, until recently). 
There is a proof of convergence in the LP'- sense in [S] . There is also a recent proof of the point- 
wise convergence of the limit within the context of uniform convergence in ergodic theorems 
in [12]. From the point of view of our present paper, it is these functions which lie at its 
heart. It is interesting to note that the way in which the isometric embedding connection 
between diffraction and dynamics is defined, the eigenfunctions are nowhere in sight. It is 
only in the L^-completion that the eigenfunctions appear, and even then the way in which 
they map (by the Bombieri- Taylor formula) is nothing like the original defining map, and has 
to be proved. 

Now we come to our present paper. Of course the underlying concern of much of diffraction 
theory is that the inverse problem (resolving structure from the diffraction) has no unique 
solution in general. The problem is that the embedding of diffraction into dynamics is not 
surjective. This problem is exacerbated in the aperiodic case. Our setting is an ergodic 
dynamical system of uniformly discrete point sets which are pure point diffractive almost 
surely. The main result of our paper is to show rather precisely the significant role that 
extinctions (places in the dynamical spectrum where there are no Bragg peaks) play in this 
ambiguity. 

In ^24J D. Mermin made the remarkable suggestion that the second and third correlations 
should always determine the structure. The argument made there is quite simple - even 
trivial - except that it again deals with quantities facing the same problem originally pointed 
out in [13] ; they are not convergent (see also Q . The author was well aware of the difficulties 
of his argument: he writes "Even granting that I have begged the question of when the density 
has a Fourier transform, and when the auto- correlation functions exist, this informal Fourier 
space argument that the identity of all second and third order correlations implies the identity 
of all higher order correlations is disarmingly trivial. I would very much like to learn of a 
comparably simple informal argument or an instructional counterexample in position space. " 
Our paper is a response to this question. If there are no extinctions, then the second and 
third correlations do suffice. If the extinctions are not too bad, we can at least get away 
with knowing only finitely many types of point correlations. As further discussed in Section 
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[5] there are also recent results which show that in fact there are situations where one really 
does need higher moments than just the second and third to resolve aperiodic point sets (in 
fact, even model sets), [TO]. 

Apart from its mathematical interest and the potential directions for further development, 
the results of our present paper seem to be physically relevant. The detailed atomic structures 
of quasicrystals are basically unknown in spite of over 20 years of work by theoreticians and 
experimentalists. Model sets are one of the primary modelling devices in the subject and 
diffraction is a fundamental tool. It is relevant to know the controlling influences on diffraction 
and to know how close diffraction, particularly in model sets, can come to determining the 
underlying structure. 

1.3. Hulls. The basic objects of interest in this paper are pairs (X, /u), where X is a set of 
r-uniformly discrete point sets A of some real space M'^ for some r > and is a probability 
measure on X. The assumptions that we need to make on {X,fj,) are listed in A(i), A(ii), 
A(iii) below. 

Throughout the paper Cr, R> 0, denotes the open cube {-R/2,R/2Y CW''. 

There is a uniform topology (called the local topology) on the set D^iW^) of all r-uniformly 
discrete subsets of M"'. The uniformity is generated from the collection of all sets (entourages) 
of the form U{K,e), K CM."^ being compact, and e > 0, where 

(1) U{K, e) = {{A, A') G VriR"^) x P^(M'^) : An K C A' + and A' n K C A + C,} . 

Thus A and A' are 'close' if on some ('large') compact K and for some ('small') e > 0, the 
points of A that are within K also lie in the e-cubical neighbourhoods of the points of A' , and 
vice- versa. It is relatively easy to see that VriM!^) is compact in this topology [29] and that 
the translation action T: 

T^{A) ■.= x + AeX, for ah x eR"^ 

of M'^ on it is continuous. An alternative description of this topology using the functions Nf 
below can be found in [5j. 
We assume 

A(i) X is a closed translation invariant subset of T)r{R'^)', 
A(ii) fi is an ergodic probability Borel measure on X. 

A(i) obviously implies that X is compact and, together with A(ii), the pair (X, /u) along 
with the group action of M'^ by translation is a dynamical system, both in the topological and 
measure theoretic senses. The assumption of ergodicity is that n is a translation invariant 
measure and X cannot be decomposed into two measurable invariant subsets which each have 
positive measure. We can, if we wish interpret as a measure on VriW^) whose support lies 
inside X. This makes it clear that ^ is the actual relevant piece of data. The space X is only 
noted for convenience. 

The basic open neighbourhoods of A defined by the uniformity on X are of the form 

U{K,e)[A] := {A' (£ X : {A, A') G U{K,e)}. 

These consist of the point sets A' that are sufficiently close to making the same pattern as 
A within the compact set K C W^. The interpretation of iJ.{U {K , e)[A]) is that it is the 
probability that a random element A' of X will lie in U{K,e)[A]. The measure fj, thus gives 
the information about what patterns are possible and what their probabilities of occurrence 
are. Its support specifies which subsets of VriM!^) are relevant. The ergodicity says that. 
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when viewed from the origin, the translations of any element A from the support of /x will, 
almost surely, faithfully represent all possible local patterns with the correct frequencies. 

We take the attitude that this is all we can hope to know about our physical system, and 
thus our objective is to determine fj, from other, physically observable, data. In our case this 
other data will consist of various correlations of the system {X, fi) . 

1.4. Diffraction and pure pointedness. Let 5'(M'^) denote the Schwartz space of all complex- 
valued infinitely many times differentiable rapidly decreasing functions on the real space M*^. 

For each n = 1, 2, . . . the n + 1-point correlation of yl G ^^{W^) is the measure 7^^^^^^ 
on M'^ X • • • X M'' (n-factors) defined by 

7^a'''Hf) = lim E ^(-^ + yi,...,-x + yn) 

R-*oo vol Cr , 

x,yi,...,yn€AnCR 

for all F G S{R'^ x • • • x M'^), if this hmit exists. 

Theorem 1.1. [51 HI [8] Let{X,fi) satisfy A{i) and A{u). Then 

(i) for fi almost every A £ X, all of the n-point correlations 7^"'* exist. Furthermore, they 
are almost surely independent of the point-set A chosen in X; 

(ii) the 2-point correlation is almost surely Fourier transformable. 

The common n-point correlations are denoted simply as 7^"\ and even more simply as 7 

(2) ^ 
for n = 2. The Fourier transform of 7^ is almost surely the Fourier transform 7 of 7. 

Definition 1. 7 is the diff"raction of (X, ^u). 

Starting with the work of Hof [13] the rigorous mathematical study of diffraction for ape- 
riodic order has attracted quite some attention in recent years. We refer to [16^ \T9\ [20] for 
recent surveys. 

A basic idea is that the correlations are, in principle, quantities that can be physically mea- 
sured. Certainly measurement of the diffraction is standard, and hence its Fourier transform, 
the 2-point correlation, may be considered as known. There are reports of inference of higher 
correlations through fluctuation microscopy |:33j , though it is not clear that these correlations 
go beyond pair-pair correlations arising from squaring the 2-point correlation. 

The space (X, //) of square integrable functions on X gets an M*^ action through trans- 
lation of functions: (Ttf){x) = f{—t + x). A simple consequence of the translation invariance 
of n is that this action is unitary for the basic inner product (•, •) defined by 

(/, 9) = fgdfi . 
Jx 

As usual / € L^(X, ^) with / 7^ is an eigenfunction of T (to the eigenvalue k E K'^) if 
Ttf = exp(— 27riA; • t)f for all t G M*^. The fact that is assumed ergodic implies that the 
multiplicity of each eigenvalue is one. The dynamical system {X, ^) is said to be pure point 
if L'^{X,fj,) has a Hilbert basis of eigenfunctions. 

A key point is a theorem that relates L'^(]R'^,7) and L'^{X,jx) and then relates the two 
concepts of pure pointedness. In order to discuss this further we need some more notation. 
For each / G S{W^) \ei Nf : X be defined by 

Nf{A) = Y,f{x)- 
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Lemma 1.2. The algebra generated by the Nf, f € 5'(R ), is dense in the algebra of contin- 
uous functions on X equiped with the supremum norm. In particular, it is dense in L'^{X,fi). 

Proof. (See [2] for a similar argument.) Obviously, the algebra in question separates points, 
is closed under taking complex conjugates and to each A £ X, there exists an / € S'(]R'^) with 
Nf{A) ^ 0. Thus, the first statement follows by the Stone- Weierstrass theorem (see [30j for 
the version used here). The last statement is then clear. □ 

Define an action U of M'^ on L^(M'^,7) by 

(C/J)(x) = e"2-*-7(x) 

for ah t,x G M'^, / e L^{R'^,j). 

Theorem 1.3. Let {X,fi) satisfy A(i) and A(ii). Then the following holds. 

(i) The set {f : f e SiW^)} is dense in L'^iR'^,^) and there is a unique isometric embed- 
ding 

e:L\R\^)^L\X,p.) 

with 9{f) = Nf for all f E S'(M'^). This embedding intertwines U and T. 

(ii) ^ is a pure point measure if and only if {X, /u) is pure point. 

(iii) For A; G M'^ the equation 

holds, where the limit is meant in the sense. Moreover, "y{{k}) 7^ i/ and only if 
0{lk) 7^ 0. In this case, 9{lk) is an eigenfunction of {X,iJ,) for the eigenvalue k. 

Remark 1. (a) The theorem may be found in the form stated here in [8]. It has a long history 
that includes [HI |T7l [121 [21 [18] . fact, (i) and (iii) are shown in [8], see [18] for extensions of 
(iii) as well. The statement (ii) is proven in various levels of generality in [11 1 117 1 fl^ l [2l [22]. 
We will give an independent proof below in Corollary 12.21 

(b) Under 6 eigenfunctions go to eigenfunctions. However, the formula of part (i) is not 
applicable in (iii): the function which takes the value 1 at and zero everywhere else, 
is not even remotely in 5'(M'^). The limit stated here appears only after approximation by 
functions from S'(IR'^). 

(c) As 6 is an isometry, j{{k}) ^ is obviously equivalent to 0(1^) 7^ 0. 

The functions ^(Ifc) for k with j{{k}) 7^ appearing in (iii) of the previous theorem will 
play a crucial role in our considerations. We define 

fk := 0{lk) whenever ji{k}) / 0. 

Definition 2. Let 

£ := {keM.'^ : k is an eigenvalue of {X, fi) } 
S := {keR'^ : 7({A:}) ^ 0}. 
£ is the dynamical spectrum of {X, fi) and S is its Bragg spectrum. 

We note that the Bragg spectrum is sometimes also known as Fourier-Bohr spectrum. For 
convenience of notation, we will write 'j{k) instead of 'j{{k}) in the remaining part of the 
paper. 
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Let us briefly discuss the structure of 8 and S and their relationship. It is well known 
that £" is a subgroup of W^. In fact, the product of eigenfunctions is again an eigenfunction 
to the sum of the respective eigenvalues and the complex conjugate of an eigenfunction is an 
eigenfunction to the inverse of the corresponding eigenvalue. 

As shown in (iii) of Thm. 11.31 any Bragg peak k comes with a canoncial eigenfunction 
/fc = ^(Ifc). In particular, we have the inclusion 

SdE. 

In general, this inclusion is strict even in the pure point case. The limit formula for in 
Thm. 11.31 shows that 

(2) f-k = Tk 
and that 

/o/O. 

Thus, the statement on the eigenfunctions in (iii) of Thm. [L3l shows that S satisfies 

G 5, and 5 = -S. 

It is a fundamental fact that the canonical eigenfunction of the Bragg peak k can be 
related to the intensity of the Bragg peak. More precisely, note that due to the ergodicity 
the modulus of any eigenfunction is constant /^-almost everywhere. Thus, 

(3) \fk\ = {fk, A>i2^(x,^) = (Ifc, lfc)L2^(Rd,=^) = likY^^ , 

where the first equality holds //-almost everywhere. 

Our final assumption is: 
A (iii) 7 is a pure point measure. 

Definition 3. A pair {X, ji) satisfying axioms A(i),(ii),(iii) is called a pure point ergodic 
uniformly discrete point process. 

In this pure point case, S generates £" as a group, as shown in [2]- We will give an 
independent proof based on Thm. 11.31 in Corollary 12.21 below. 

Definition 4. The set X := £\S is called the set of extinctions of {X,fi). 

1.5. Correlations and moments. Let {X,fi) satisfy A(i) and A(ii). For n = 1, 2, . . . , the 
nth-moment of fi is the measure /i„ defined on M'^ x • • • x R'^ (n-factors) by 

fin{hi, . . . ,hn) = Nhj^ . . . Nh„ djj. 

Jx 

for all hi,...,hn E S'(M'^). It is clear that these moment measures are invariant under 
simultaneous translation of all the variables. There is a standard procedure of eliminating this 
translation invariance resulting in the reduced moments fx^^'^ which are in one less variable. 
This works as follows: let g,hi, . . . , hn-i € S{W^) and let G be the function on {W^Y whose 
value on (x,?/i, . . . ,yn-i) is g{x){T^hi){yi) ■ ■ ■ {T^hn-i){yn) ■ Then 

fin{G) = fin{g{Txhi) . . . (T^hn-i)) = / g{x)dx fi'^'^ {hi, . . . , hn-i) , 

and this equation defines the reduced moments. 
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Most importantly for our purposes, these reduced moments are also the correlations - see 
[7], Sec. 12.2 and [8], Sec. 7u- More precisely, the following holds. 

Theorem 1.4. Let {X,fi) satisfy A(i) and A(ii). 

(i) For each m G N, /i zs uniquely determined by its moments n > m. 

(ii) 7("^=M^^^ri = 2,3,.... 

(iii) For n > 2, is uniquely determined by /i^'^'^. 

(iv) The measure is uniquely determined by n = 2, 3, . . . . 

Proof, (i) For m = 1 the statement follows immediately from Lemma 11.21 Now, it suffices 
to show that the /U„, n > m, determine fim- By Lemma ll.2| again, the constant function 1 
can be approximated by elements of the algebra generated by the Nfi, h G S'(R'^). Thus, a 
product 

N,,,---Nh^=Nh,...Nj,^.l 

can be approximated by linear combinations of products of more than m functions in S{M.'^). 
Thus, fin, n > m, determine fim- 

The proof of (ii) can be found in [7J, Prop. 12. 2. V. The proof of (iii) can be found in [7], 
Sec. 10.4. Finally, (iv) is a direct consequence of (i), (ii) and (iii). □ 

Remark 2. The proof of (i) in the previous theorem does not require ergodicity. It only uses 
that the functions Nj, f £ S{W^), are bounded and continuous on the compact X. 

The point of the previous theorem is that rather than correlations, we may instead look at 
corresponding moments. Our question becomes that of asking how many moment measures 
are required to pin down ii uniquely. 

2. ElGENFUNCTIONS AND CYCLES 

2.1. The cycle function of (X, ^u). Let {X,fi) be a pure point ergodic uniformly discrete 
point process. Thus, A(i),(ii),(iii) are valid. 

The elements of L'^{R'^,^) are all the sums 

Xfclfc, where ^ \xk\'^^{k) < oo . 
fees kes 

As described in (i) and (iii) of Thm. 11.31 the map 9 exhibits very different behaviour on 
functions h € S'(M'^) and on functions 1^. This leads to two very different ways in which to 
write ^(L^(M'^, 7)). More precisely, both the linear span of the set of 1^, k € S, and the set 
h, h £ 5(M'^), are dense in L^(M'^,7). As 9 is an isometry, this gives 

9{L^{R'^, 7)) = {Nh : h£ S{R'^)} = linear span{/fc : k £ S} . 

Our next aim is to obtain similar statements for products. This requires some care as we will 
have to deal with products of infinite sums. The corresponding details are given in the next 
two lemmas. 

We will use repeatedly the elementary fact that {gmh} converges to gh in whenever 
{gm} is a sequence converging to g in and /i is a bounded function. 

^The reduced moments are also connected directly to Palm measures, a direction that is more fully explored 
in [miH]. 
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Lemma 2.1. Let n € N and /ii, . . . , /i„ € S'(M°') be given. Then, 

kies kn&s 

where the sums exist in and are taken one after the other. In particular, 

J Nh^... Nh^dn = E • • • E ^l(^l) • • • fT'n{kn)lJ'{fki ■ ■ ■ fkj- 
kiGS kn&S 

Proof. By (i) of Thm. 11.31 we have 

Nh, = J2 hj{k,)fk, 

for each j. Therefore, 

Nh,... Nh„ = I E ^lih) fk, \Nh,...Nh^=J2 hi{ki)fk,Nh, . . . iV/.„ 

\fciG5 / fciG5 

and the first statement follows by induction. As /i is a finite measure, the last statement then 
follows easily. □ 

As a corollary of this lemma we obtain a new proof of the following known fact. 

Corollary 2.2. T has pure point spectrum, i.e. there exists a basis of L'^{X,fi) consisting of 
eigenf unctions. Moreover, any eigenfunction is a finite product of functions fk, k ^ S, and 
any eigenvalue is a sum of k G S. 



Proof. By the previous lemma, any function of the form N^^ . . . Nh„ with hj € ^(IR ) can 
be approximated by linear combinations of products of the form /^^ . . . f^^, kj G S. Lemma 
11.21 then gives that functions of the form fk^ ■ ■ ■ fk„, n G N, /cj G 5, j = 1, . . . , n are total 
in {X, fi) . As each function of the form fk^ ■ ■ ■ fk„ is an eigenfunction to the eigenvalue 
ki + ■ ■ ■ + kn the statement follows. □ 

Lemma 2.3. Let n G N and ki, . . . ,kn G 5 be given. For each j = 1, . . . ,n, let {hj^^} be a 
sequence in S{W^) whose Fourier transforms converge to Ik- in L^{W^,j). Then, 

/fci • • • fk„ = lim lim ... lim (mj) . . . A , 

mi^oom2^oo m„— ♦oo "j ™n 

where the limits are taken in L^. In particular. 



Kfki ■ ■ ■ fkj = Jim_ Jim_ . . . _lim_ n{N^(mi) ■ ■ ■ 



Proof The functions N^, and the functions fk, , = 1, . . . , n, G N are bounded. Thus, 

the convergence of the /i^™"^ easily yields convergence of the products. (Note that the limits are 
taken one after the other). As /u is a finite measure the last statement then follows easily. □ 

We will now introduce a crucial object in our studies, namely the cycle function a. Notice 
that, by almost sure constancy of the modulus of the functions fk, 

l/fci ■ ■ • /fc„ P = / /fci • • • fk„fkn ■ ■ ■ fkidfl 

Jx 

= {fkiJk^)---{fk„,fkJ = 7(/ci)...7(A;n). 
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If fci + • • • + A;„ = then fk^--- fk„ is an eigenvector for 0, and hence is a multiple of the 
constant function Ix- Thus, in this case, 

(4) fk,-.. fk^ = aih, k^mh)^/^ . . . 7(A:„)V2 
for some 

a{ki,...,kn) G U{1). 

Here, U{1) is the unit circle i.e. the set of all complex numbers of modulus one. For any 
ki, . . . ,kn € S, we therefore obtain 

(5) fJ-{fkr---fk,J = / /fci • • -Anlx^M = (/fci • • -/fcnjlx) 

Jx 

('a(fci,...,fc„)7(fci)V2...^(A:„)i/2 a k, + . . . + k^ = 0; 

\o ifki + --- + kny^O, 

since eigenf unctions for different eigenvalues are orthogonal. 

The next two results basically say that knowledge of the cycle function a determines the 
moments and vice versa. 

Proposition 2.4. Let n G N 6e given. Then, the n-th moment of fi is uniquely determined 
by 7 and the quantities 

a{ki , . . . ,kn) for ki, . . . ,kn € 5 with ki + ■ ■ ■ + kn = ■ 

Proof. This follows directly from Lemma |2. II and ([5|). □ 

Proposition 2.5. Let n G N 6e given. Then, the values a{ki, . . . , /c„) for ki, . . . ,kn G 5 with 
ki + ■ ■ ■ + kn = are uniquely determined by 7 and the n-th moment of 

Proof. This follows from Lemma 12.31 and ([5]) . □ 

It is convenient to introduce the Cayley graph ^ of £^ with respect to the set of generators 
S. Its vertices are the points of £ and its edges are the pairs {k, 1} of vertices whose differences 
A; — Hie in S. Since S = —S, we may treat the edges as undirected. Any k = (/ci, . . . kn) G 5" 
with ki + . . . + kn = leads to a cycle {0, ki, ki + k2, ■ ■ ■ , ki + ■ ■ ■ + kn-i, fci + • • • + A;„ = 0} in 
Q. Thus the function a described above can be thought of as a function of the set Z of cycles 
of G- We shall call it the cycle function of {X,fj,). 

2.2. Properties of the cycle function. Given k, 1 G their concatenation 

kl := (fci, . . . , kji, li, ■ ■ ■ , Ip) , 

is obviously also in Z. 

Proposition 2.6. The cycle function a has the following properties: 

(i) for all k,leZ, a(k)a(l) = a(kl); 

(ii) a(0) = 1; 

(iii) for all \<. Z, a(k, — k) = 1; 

(iv) a{ki, . . . , kn) is independent of the order of the elements ki,...,kn making up the 
cycle; 

(v) given any cycle \<. (z Z, then any pair {k, —k} where k ^ S, and also can be inserted 
into or deleted from the symbols of k without affecting the value of the cycle function 
a. 
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Proof, (i) follows from (g]). As for (ii), note that /o = 0(0)7(0)1/2 by dH). Since /o = 6i(l{o}) > 

from Tlim.[L3](iii), 7(0) > 0, and a(0) e ^7(1), we see that a(0) = 1. 
This proves (ii). Part (iii) follows from 

a{k, -k)^{k) = fkf-k = fkfk = l/fcP = l{k). 
Items (iv) and (v) are trivial consequences of (jl]) and parts (i),(ii), and (iii). □ 

Let Zn ■= {{ki, . . . , kn) G Z}, Zq := {0}, and Z^o = \J^=o Zn- We introduce an equivalence 
relation on Z^o by transitive extension of the two rules: 

• k ~ 1 if 1 is a permutation of the symbols of k 

• k ~ 1 if 1 can be obtained from k by inserting or removing pairs {k, —k}, k £ S, 01 
by inserting or removing 0. Let 

Z := Zoo/ ~ • 

It is easy to see that k ~ 1 , k' 1' =^ kl ~ k'l' , so multiplication descends from Zoo to Z. 
Indeed Z is an abelian group under this multiplication, with 0~ as the identity element. 

Of the various {ki, . . . , kn) that can represent a given element k (z Z there is (at least) one 
of minimal length n. This minimal length is denoted by len(K). Define 

Zn := {k, Z : len(K) < n} . 

We shall also write len(k) = len(K) when k"" = k and call it the reduced length of k. 

Z = \JZn, 

n=0 

ZnZp C Zn+p for all n, p. 
Evidently the cycle function a determines a homomorphism, a, 

d : Z — > U{1) 

with a(K) = a{ki, . . . , kn) if k = (A;i, . . . , kn)^ ■ It is clear from this that d is known on Zg by 
its values on the sets Zn for n < q if 

Zq = Zn Zp . 

n+p=q,0<n<q 

2.3. Main results. 

Theorem 2.7. Let {X, fi) be a pure point uniformly discrete ergodic point process and suppose 
that 

S + ■■■ + S = £ . 

n 

Then Z is generated, as a group, by ^2n+i ■ 

Proof: Let k = (ki, . . . , k^) be any cycle where N > 2n+l. By assumption ki-\-- ■ ■+kn+i G S 
can be written in the form + • ■ ■ + where the li E S. Then with j := (ki, . . . , kn+i) and 

1 = {li,. . . ,ln) 

k~((j)(-l))(l((-j)k)) 

which writes it as the product of (j)(— 1) and 1((— j)k). These have reduced lengths at most 
2n + 1 and A'" — 1 respectively, and this shows that Zjy C Z2n+iZiy-i. The proof finishes by 
induction. □ 
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Theorem 2.8. Let {X, fj,) be a pure point uniformly discrete ergodic point process. If the 
dynamical spectrum is finitely generated by the diffraction spectrum, then fj, is uniquely deter- 
mined by a finite number of its moments ^rn- More precisely, if 

S + -.. + S = £ 
^ V ' 

n 

then n is uniquely determined by its moments fim, m = 2, . . . , 2n + 1. 

Proof. By (ii) of Thm. 11.41 the second moment determines the autocorrelation 7. By Propo- 
sition 12.51 the moments /imj = 2,...,2n-|-l, then determine the function a on for 
m = 2, . . . ,2n + 1. The previous theorem gives that the function a is then completely deter- 
mined. Thus, by Proposition 12.41 all moments fim, m > 2, are then determined. Now, the 
theorem follows by (i) of Thm. 11.41 □ 

Corollary 2.9. If there are no extinctions, so S = £, then // is determined by its second and 
third correlations. 

Proof. This follows immediately from the previous theorem and parts (ii) and (iii) of Thm. fL^ 

□ 



3. Model sets 

The theory of model sets is a good place to find examples of the types of point processes 
{X, |u) that we have been discussing. We start with a cut and project scheme {W^, H, L) with 
corresponding 'torus' T := iW^ x H)/L. We suppose that H \s a. complete metric space. The 
canonical mapping from the projected image L of L in to its projected image in H is 
denoted by (•)*. We suppose also that we have a regular model set A = yl(VF) := {x G L : 
X* G W} given by some regular closed set W d H (i.e. W = W° and W has boundary of 
measure 0). In this case, the hull X = X{A) of A is uniquely ergodic, all elements of X share 
a common autocorrelation 7, and the diffraction 7 is pure point [31]. See [251 13] for basic 
material on model sets. 

We wish to consider the situation regarding the diffraction and extinctions. Let (M'^, H, L°) 
be the corresponding dual cut and project scheme, where L° is the dual group of T, [2S], Sec. 5. 
The natural projected image ofZ° in M'^ ~M'=' is denoted by £. We shall also use the ★-notation 
for the dual cut and project scheme. Since the projection of L° into H has dense image, £* 
is dense in H. 

The diffraction of A is known [14', '31] to be given by 

(6) m = \hvi-k*)\\ 

for all k (z £. Here, the measure of the torus T is normalized to be one. 

The set of extinctions X is then the set of k € £ for which liy(— A;*) vanishes, and the 
Bragg peaks make up the set S = £\X. We know that S generates £ as a, group and G 5. 

Lemma 3.1. Suppose that X* has no interior. Then S + S = £. 

Proof. We need to prove that every element of X is the sum of two elements of S. Let K be 
the kernel of ( )*. Prom the form of the diffraction Q, both S and X consist of unions of 
full cosets of K. Thus it suffices to show that every element of X* is the sum of two elements 
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of S*. Let z G X* and suppose that z is not so expressible, i.e., {z — 5*) n 5* = 0. Then 
z-S*CX*,soS*Cz-^. From this 

H = £* = S*\JX* = S^u7^c{z-'X*)u'X*, 

which is impossible by the Baire category theorem (since X* has no interior). □ 

Remark 3. One should note that even in the same cut and project scheme, different windows 
can give rise to the same diffraction in ([6]), even though the windows are not translationally 
equivalent. An example of this, that derives from the covariogram problem, is discussed, 
along with references to covariogram literature, in [1]. This shows that the second moment 
alone is insufficient even to distinguish model sets from the same cut and project scheme. 

4. Origins of the problem in questions of symmetry 

The results in this paper have a number of points of contact with the ideas of D. Mermin 
[23] and subsequent works of R. Lifshitz, D .A. Rabson, and B. N. Fisher [23l|26]. The start- 
ing point was a puzzle which arose almost immediately after the discovery of quasicrystals. 
The usual notions of symmetry from crystallography are not adequate in the theory of qua- 
sicrystals. First of all, translational symmetry is drastically diminished, often to the point of 
being non-existent. Second, even the finite symmetries are somewhat nebulous. Tcosahedrally 
symmetric' quasicrystals, for example, have perfectly isosahedrally symmetric diffraction pat- 
terns, but they need not be literally icosahedrally symmetric in the sense that the structure 
is mapped precisely onto itself by the point symmetries of the icosahedral group. In this 
section we will discuss this question and how it can be resolved through the use of dynamical 
systems. 

Mermin's solution to the symmetry problem was to realize it on the Fourier side of the 
picture rather than the physical side. Briefly, the idea was to express the density distribution 
p of the quasicrystal as a superposition of plane waves from its translation module (what we 
have called £ above): 

(7) p(r) = ^/?(fc)e2-^-, 

and then to remark that two densities p and p' based on the same module of wave vectors 
are physically indistinguishable if their correlations of all orders are identical, something that 
should happen if 

(8) p{ki)...p{kn)=p'{ki)...p'{kn) 

for all A;i, . . . fen G £^ with fei + • • • + A:„ = 0. Thus symmetry becomes symmetry in the sense 
of indistinguishable correlations. 

Now yO and p' will be indistinguishable if for some x • £ — ^ M/Z we have 

(9) p'{k) = e2^*^(^)p(A;) . 

and a symmetry g of £ would appear as a symmetry of the physical system if, for all k £, 

p{g{k)) = e^^'^^^^W) 

for some suitable Xg '■ ^ ^ M/Z. 

Applying ([8]) at n = 3 in ([9]) already gives 

xih + k2) = xih) + xih) 
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for all ki,k2 € £, which already provides all the information derivable from ([8]) for all other 
n. It was from this that Mermin concluded that 2- and 3-point correlations should determine 
everything. 

The main difficulty with this approach, and this is already made clear in [23] , is to give any 
mathematical meaning to the expressions ([7]), except as tempered distributions. This issue 
is discussed in some detail by A. Hof in [Hj. The measure p representing the distribution 
of density of the quasicrystal, say p = ^^eA^^^ ™ general Fourier transformable 

as a measure. If instead one treats it as a distribution then it is difficult to say what its 
Fourier transform looks like and particularly whether or not it is composed of a countable 
sum of weighted deltas at some (usually dense) subset of M.'^. Certainly the formation of 
moments in this language would be a formidable task. The upshot of Hof 's study of diffraction 
was his approach to it using the Fourier transform of the autocorrelation (which is Fourier 
transformable as a measure), and this has been the basis of most subsequent mathematical 
work on diffraction. Thus |24j is better seen as a formal vision of how things should work out 
rather than a rigorous exposition of how they actually do, and in that sense it prescient. 

In our approach the sums ^ appear in a form that is tamed by averaging, namely the 
expressions appearing in Thm. I1.3l (iii). Such an expression is non-zero only if A; € f is 
the position of a Bragg peak, i.e. there is no extinction there. In [24j it is claimed that 
extinctions can only occur due to symmetries, but the examples from model sets show that 
extinctions come from the Fourier transform of the window function and do not seem to be 
related specifically to symmetries. Furthermore the extinctions seem to be potential, but not 
absolute, obstructions to the correctness of the assertion about 2- and 3-point correlations, 
see ^ below. 

In hindsight one can see that the symmetry question amounts to symmetries of the cor- 
relations, hence of the moments, and ultimately of the measure p itself. An isometry g of 
gives rise to a mapping A i— > gA. This g is a symmetry of the point process if p is 
(7-invariant. If X is the support of p, then this also entails that g{X) C X. This is the idea 
of symmetry put forward by Radin in [27^ [28] . 

This symmetry can also be expressed by means of groupoids. In fact, there is a canonical 
groupoid structure associated to the dynamical system (X, M'^). Associated with it is the 
groupoid of the transversal (see [5] and references therein). The point groups of symmetries 
of the system then act as isomorphisms of this transversal groupoid. The measure p appears 
by giving a trace on the corresponding C*-algebra. In this picture, invariance of the system 
under a symmetry means invariance of the trace under this symmetry. This invariance of 
the trace then amounts to invariance of the measure p, which in turn implies invariance of 
diffraction and higher moments under the symmetry. 

5. Final comments 

There are a number of obvious questions that arise from this work. Foremost is the question 
of whether every cycle on Z arises from some sort of dynamical system of density distributions 
on W^. In [21] we have taken a first look at this. This requires changing the setting from the 
study of uniformly discrete subsets of M'^ to more general distributions of density arising by 
formalizing the main ideas of this current paper. 

There is also the interesting question of the real role of extinctions. Although we have seen 
that extinctions are an obstruction to our method of reconstructing the measure p from its 
moments, we have not shown that this obstruction is one of principle rather than an artifact 
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of our approach. Inspired by the present paper this topic has been taken up recently: In [9] 
it is shown that for model sets with real internal spaces, whether or not there are extinctions, 
the 2- and 3-point correlations determine the model set within the class of all model sets. 
On the other hand, as soon as the internal space is not purely real model sets exist in which 
the extinctions force higher correlations to be used. The paper |10j offers examples of multi- 
atomic model sets in which the scattering strengths of the different point types are differently 
weighted and for which even the 2, 3, 4, 5-point correlations are insufficient to resolve them. 

Another question is the curious role of odd numbers in our main result. Are there general 
conditions under which even numbered higher correlations play the defining role? 

Finally it would be interesting to develop practical methods for actually constructing in- 
creasingly accurate atomic approximations to the densities using the increasing information 
from the 2-, 3-, ... correlations. 
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